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ABSTRACT 

We present estimates for the size and the logarithmic slope of the disk tem¬ 
perature profile of the lensed quasar Q2237-I-0305, independent of the component 
velocities. These estimates are based on six epochs of multi-wavelength narrow- 
band images from the Nordic Optical Telescope. For each pair of lensed im¬ 
ages and each photometric band, we determine the microlensing amplitude and 
chromaticity using pre-existing mid-IR photometry to define the baseline for no 
microlensing magniheation. A statistical comparison of the combined microlens¬ 
ing data (6 epochs x 5 narrow bands x 6 image pairs) with simulations based 
on microlensing magniheation maps gives Bayesian estimates for the half-light 
radius of R 1/2 = (M)/0.3 Mq light-days, and p = 0.95 ± 0.33 for the 

exponent of the logarithmic temperature prohle T oc This size estimate is 

in good agreement with most recent studies. Other works based on the study of 
single microlensing events predict smaller sizes, but could be statistically biased 
by focusing on high-magniheation events. 
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Subject headings: accretion, accretion disks — gravitational leasing: micro — 
quasars: individual (Q2237+0305) 


1. Introduction 

The basic model for describing the inner regions of quasars is the thin disk model 
(Shakura & Sunyaev 1973; Novikov & Thorne 1973) which predicts the size of the accretion 
disk and the radial dependence of its surface temperature. Gravitational microlensing (Chang 
& Refsdal 1979, 1984; see also Kochanek 2004 and Wambsganss 2006) is the main tool used to 
estimate both parameters, either from time variability or through the wavelength dependence 
of the microlensing magnihcation. Microlensing studies (see e.g. Pooley et ah 2007; Morgan 
et ah 2010; Sluse et ah 2011, Blackburne et ah 2011, 2014, 2015; Mosquera & Kochanek 
2011; Jimenez-Vicente et ah 2012, 2014; Hainline et ah 2013; Mosquera et ah 2013; MacLeod 
et ah 2015) have found that the mean sizes of quasar accretion disks are roughly a factor 
of 2-3 greater than the predictions of the standard thin disk model. These differences are 
too large to be explained by contamination from the broad emission lines and the pseudo¬ 
continuum contributions, or scattering on scales larger than the accretion disk (Dai et al. 

2010, Morgan et al. 2010). Recent measurements of wavelength-dependent continuum lags in 
two local active galactic nuclei (AGNs) are consistent with the microlensing results (Shappee 
et al. 2014, Edelson et al. 2015, Fausnaugh et al. 2015). 

The study of the disk temperature prohle is more complicated because multi-wavelength 
observations are needed to detect chromatic microlensing (Anguita et al. 2008; Poindexter 
et al. 2008; Bate et al. 2008; Eigenbrod et al. 2008; Floyd et al. 2009; Blackburne et al. 

2011, 2015; Mediavilla et al. 2011a; Mosquera et al. 2011; Munoz et al. 2011; Motta et 
al. 2012, Rojas et al. 2014). In this case, it also becomes more important to separate the 
contributions from strong emission lines. Such a separation can be done most cleanly using 
narrow band photometry or spectroscopy. 

Many of the studies of the source size in Q 2237-1-0305 are based on the htting of its light 
curves with tracks extracted from simulated microlensing magnihcation maps (Kochanek 
2004; Poindexter & Kochanek 2010b; Mosquera et al. 2013). As an alternative to light-curve 
htting, we will study the size and temperature prohle of Q 2237-1-0305 using several epochs 
of multi-wavelength narrowband observations. This allows us to remove the inhuence of the 
broad emission lines on the amplitude of the continuum microlensing (Mosquera et al. 2009). 
We will follow a procedure similar to that used in SBS 0909-1-532 (Mediavilla et al. 2011a) 
and HE 1104-1804 (Munoz et al. 2011), but in the case of Q 2237-1-0305 we have signihcantly 
better statistics with six epochs (four will be considered independent) and four quasar images. 
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This method relies on changes in the microlensing amplitnde and chromaticity bnt not on 
its dependence with time and hence no velocity estimates are necessary. However, a baseline 
for the intrinsic flnx ratios in the absence of microlensing is needed, and we will dehne these 
by the mid-IR flnx ratios from Minezaki et al. (2009, see also Agol et ah 2000). In Section 
2 we describe the observations and data. Section 3 is devoted to the statistical analysis and 
in Section 4 we discnss the resnlts. 


2. Observations and Data Analysis 


Q 2237+0305 (Hnchra et al. 1985) was observed with the 2.56 m Nordic Optical Tele¬ 
scope (NOT) located at the Roqne de los Mnchachos Observatory, La Palma (Spain). We 
nsed the 2048x2048 ALSFOC detector, with a spatial scale of 0.188 arcsec/pixel. We ob¬ 
tained a total of five epochs in 2006 and 2007 which we combine with an earlier epoch from 
Mosqnera et al. (2009) taken in 2003. A set of seven narrow Liters plus the wide Bessel I 
Liter were used, covering the wavelength range 3510-8130 A. Table [1] provides a log of our 
observations. For the Lrst epoch, a V-band image was also taken. In the second epoch, the 
wavelength coverage was poorer, because we were unable to observe with all the Liters due 
to bad weather conditions. 


The data were reduced using standard irafQ procedures, and pont-spread function 
(PSF) photometry Ltting was used to deriv e the diLerence i n mag nitude as a function of 
wavelength between the four images. As in iMosquera et al.l (1200911 . the galaxy bulge was 
modeled with a de Vaucouleurs proLle, and the quasar images as point sources. This model 
was convolved with PSFs derived f rom stars obse r ved in each of the frarn e s, and Lt to the 
image using statistics following iMcLeod et al.l (119981) and iLehar et al.l (120001) . Even in 
the bluest Liters, the Lts were excellent due to the good seeing conditions for most of the 
epochs (076 in I-band). Among the narrower Liters that we used, only two were aLected by 
the broad emission lines of the quasar. The Stromgren-u Liter contains roughly 40% of the 


Lya e mission line, based on the SDSS quasar composite spectrum from IVanden Berk et al. 


( 2001 '). And at A = 4110 A, the wavelength range covered by the Stromgren- y Liter coincides 


with the position and width of the CIV emission line (for more details see iMosquera et al. 


20091 1. 


^IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Associ¬ 
ation of Universities for Research in Astronomy (AURA) under a cooperative agreement with the National 
Science Foundation 
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Figure [T] shows the results for the six nights and Table [T] reports the photometry for 
the new data. A clear wavelength dependence of the flux ratios (m* — mA)i=B,c,D relative 
to image A is observed for the epochs HJD 2454056 (Fig. [T]i) and HJD 2454404 (Fig. 
m). Chromatic effects are also observed relative to image B {rrii — mB)i=A,c,D on the hrst of 
these two dates. Since extinction effects would be common to all the epochs (e.g. Falco et ah 
1999 or Munoz et al. 2004) and the time delays are negligible, this depende nce can only be 


expla ined by chromatic microlensing of images A and B (for more details see iMosquera et al. 


20091 1. Since image C seems little affected by chromatic microlensing on any of those nights, 
the difference {me — mAY — {me — mAY allows us to determine the chromaticity in image 
A, with values of m^ — = —0.43 ± 0.16 and — m^ = —0.33 ± 0.16 for HJD 2454056 

and HJD 2454404 respectively. In a similar way, the chromaticity for image B is m^ — 
= —0.25 ± 0.16 for HJD 2454056. The I and b hlters were chosen to perform these 
calculations since they are the reddest and bluest hlters that are not affected by the broad 
emission lines. The other epochs do not show any signiheant dependence on wavelength. In 
particular the magnitude difference as a function of wavelength for the epoch HJD 2454001 
(Fig. Ic) is hat given our errors of order ~ 0.1 mag. This suggests a smaller diherential 
extinction than the estimates of Agol et al. (2000), although we cannot rule out fortuitous 
cancellations of color variations due to extinction by chromatic microlensing ehects. Finding 
chromatic microlensing during this period of observations is not a surp rise, since the Optica l 
Gravitational Lensing Experiment (OGLE) V-band continuum data fjWozniak et al.ll2000[) 
shows signiheant brightness variations in the three independent magnitude diherences at 
this time. This implies that microlensing is strongly ahecting the quasar images, making 
chromatic microlensing ehects more likely. 


3. Results 


The statistical procedure used to estimate the size of the accretion d i sk and the strength 
of the chromati city is a variant of tha t desc ribed in iMediavilla et al.l (l2011al) . Munoz et 
al. (2011) and Ijimenez- Vicente et al.l fl2012f ). The main diherence is that the observed 
magnitudes for the four quasar images in the hve diherent bands are used simultaneously to 
calculate the probability distributions. 

The magniheation maps were calculated using the Inverse Polygon Mapping algorithm 
(Mediavilla et al. 2006, Medi avilla et al. 20 1 1a). We have used canonical values for k 


and 7 for the four images from [ Schmidt et al. (119981) and out all the mass into equal mass 


stars. The models used bv iPoindexter fc Kochanek fcOlOal) are very similar but with a mass 
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spectrum for the stars. The 2000x2000 pixels maps have a pixel size of 0.5 light-days for a 
mean stellar mass of (M) = IMq and all linear sizes can be scaled with ({M)/. We 
considered three surface brightness prohles. In all three models, the scaling of the source size 
with wavelength is described by a power-law, r 5 (A 2 ) = (A 2 /Ai)Prs(Ai) where Aj is the rest 
wavelength corresponding to each filter. The first model we consider is a simple Gaussian, 
I{R) oc exp(—i?^/2r^). The second model, I{R) oc (exp — l) becomes the 

standard thin disk model when p = 4/3. The third model, I{R) oc (exp [{R/rsY^^] — l) ^ 
also becomes the standard thin disk model when p = 4/3. The third model corresponds to 
a thermally emitting disk with a temperature prohle T oc R~^^p and is well-dehned only for 
p > 1/2 since we are not including an inner edge. The second model is a hybrid, but by 
holding the exponent in the blackbody function fixed and only varying the scale length, the 
brightness profile is well defined for all p. We will refer to the three models as the Gaussian, 
hybrid and thin disk models. The length scales will depend on the profile and p, however, 
we expect the estimates of the half-light radii for the different profiles to agree based on the 
results of Mortonson et ah (2005). We use the bluest observed wavelength, = 4670 A as 
the reference wavelength, corresponding to Ai = 1736 A in the rest frame. When no reference 
is made to a specihc wavelength the size is for this reference wavelength, r* = (1736 A). 


The observed magnitude of image I = {A, B, C, D) is 


i) = mo(tj, f) + Pi + 6pi{tj, i) 


( 1 ) 


where mQ{tj,i) is the intrinsic magnitude of the source at time tj in hlter i, pi is the macro 
magnihcation (in mag) of image I and 6pi{t, i) is the time varying microlensing magnihcation 
(in mag) of image I in filter i. Since we have no direct information on the intrinsic variability 
of the source or the absolute magnifications, we ultimately want to work in terms of magni 


tude d ifferences between images. We assume that the mid-IR flux ratios from iMinezaki et ah 


(120091) (cf. their Table 4) correctly estimate the intrinsic flux ratios pYj, which means that 
they define all the values of = pi — pj. Following Kochanek (2004), we eliminate the 
intrinsic magnitude of the source mo{tj,i) by starting from the £t to the four individual 
images at a single epoch and wavelength 


^ - pi - 6pi{tj,i)Y ( 2 ) 

/ 

and optimizing for the best source model mo(tj,i). If we then substitute this best fit source 
model into x^(tj, i), we are left with a sum over the six possible pairs between the four images 

i) ^ [Am/(fj, i) — /Amj{tjA)f 

I .]>! 


( 3 ) 
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where 

Ami{tj, i) = - /i/ - i) (4) 

and the errors ajj(tj,i) are defined in the eqnation 7 of Kochanek (2004) bnt rednce to 
if (T/ = aj{= a) as we will assnme her^. This expression only depends 
on fij — hj, whose valnes are snpplied by the mid-IR finx ratios. Thns, the fnll Xj for a given 
epoch tj becomes the snm of Eqnation [3] over the filters, 

( 5 ) 

i 

which can be evalnated for any microlensing trial defining the valnes of 

For a given pair of parameters {rs,p), the magnification maps for the four images 
A,B,C,D are convolved to the size rs(A/Ao)^ appropriate for each wavelength and the Xj 
are calculated for N=10® randomly selected locations in each of the four magnification maps. 
The probability density function Pj{rs,p) for each epoch is computed as the sum 

N 

Fj(rs,p) oc (6) 

i=l 

over a 2D grid of values for and p. We have used a logarithmic grid in Vg such that 
ln(r*/light-days) = 0.3 x i for i = 0, • • •, 17 and a linear grid in p with p’ = 0.25 x j for 
j = 0, • • •, 9. This way, Vg spans from roughly 1 to 165 light-days and p runs from 0 to 2.25. 
For the thin disk model, which is only defined for p > 1/2, we ran extra cases at p = 0.55 
and 0.68 and then follow the remainder of the sequence. 

We have data for six different epochs between August 2003 and October 2007. Some of 
the epochs are very close on time and we have combined them into a single epoch because 
they cannot be considered independent. Epoch HJD 2453968 (Aug 21st 2006) is combined 
with HJD 2454001 (Sep 23rd 2006) and epoch HJD 2454388 (Oct 15th 2007) is combined 
with HJD 2454404 (Oct 31st 2007). The dispersion around the mean for all values of the 
combined epochs is 0.08 mags, which is the same as the measurement error we adopted for 
the merit function above. 

The joint probability distribution V{rg,p) for the four independent epochs is obtained 
by the product of the probabilities for each of the epochs 


^Although there are slightly different errors for every individual image we have chosen to use the average 
measurement error a = 0.08 for weighting all the data. 
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^{rs,p) = YlPj{r„p) (7) 

j 

The resulting joint probability distributions of the three surface brightness profiles are shown 
in Figure 2 and are compared to previous size estimates for Q 2237+0305 in Figure 3. The 
maximum likelihood estimate corresponds to V{M)/0.3Mq 

light-days andp = (0.50lo;27,1.0lo;5o, 0-85lg;25) Gaussian, hybrid and thin disk models 

respectively. Using a logarithmic (linear) prior for (p), we obtain Bayesian estimates 
for the expected values of = ( 7 . 014 ^°, 3.7^^!, 7.9l5®3®)-\/ (M)/0.3Mq light-days and p = 
(0.66 ± 0.32,1.04 ± 0.29,0.95 ± 0.33). As a test we also computed the joint probability 
distribution without merging the closely separated epochs and obtained very similar results. 


4. Discussion and Conclusions 

At Ai = 1736 A, the half-light radius estimates for the three (Gaussian, hybrid, and 
thin disk) models of R 1/2 = O.Ol®;^, 8.5l4;o)\/ {M)/O.SMq light-days are remarkably 

similar (see Figure 2), consistent with the prediction from Mortonson et ah (2005) that esti¬ 
mates of the half-light radius should be insensitive to the particular surface brightness prohle 
shape. Our R 1/2 estimations are also in agreement with previous estimates by Poindexter & 
Kochanek (2010b; R 1/2 = 5.4 ±3.2 light-days), Sluse et ah (2011; R 1/2 = light-days) 

and Mosquera et ah (2013; R 1/2 = O.Ol^g light-days), as shown in Figure[3l The latter two re¬ 
sults are also scaled to (M) = O.3M0. The Poindexter & Kochanek (2010b) result marginal¬ 
izes over the uncertainties in (M), finding (M) = 0.52Mq (O.I 2 M 0 < M < I. 94 M 0 ). 
Poindexter & Kochanek (2010ab), Mosquera et ah (2013) and Sluse et ah (2011) set the 
scales by adopting priors on the effective source velocity rather than choosing a fixed mean 
mass (M). That our results agree demonstrates that these priors are reasonable. 

Analyses of individual high magnihcation events tend to measure smaller sizes. For 
example Anguita et ah (2008) find R 1/2 = l-OlgJ light-days and Eigenbrod et ah (2008) 
hnd Ri /2 = 3.0 ± 2.0 light-days. However, studies of high magnihcation events are likely 
biased toward small quasar size estimates because it is easier to obtain high magnihcations 
with small sources (e.g. Kochanek 2004; Eigenbrod et ah 2008; Blackburne et al. 2011). It is 
also interesting to note that without considering any velocity prior (i.e. adopting a uniform 
prior), the size determinations of Eigenbrod et al. (2008) and Anguita et al. (2008) increase 
by a factor ~ 4 (see Sluse et al. 2011) but would then require cosmologically unrealistic 
peculiar velocities for the lens/source. 

In any case, most of the results derived from the optical may be reconciled near a value 





of ~ (7 ± 4) a/ (M)/0.3Mq light-days. This value is large compared to the predictions of 
the thin disk model based on the flux (~ 1 light-day, Mosquera & Kochanek 2011). It is in 
agreement with the results from other lenses, where the Morgan et al. (2010) black hole- 
mass-size correlation predicts a size of ~ 5 light-days for an estimated black hole mass of 
1.2 X IO^Mq (Assef et al. 2011). This size discrepancy is now also seen in the recent size 
estimates for two local AGN using measurements of continunm lags (Shappee et al. 2014, 
Edelson et al. 2015, Fausnangh et al. 2015). The size problem is clearly not unique to the 
microlensing method. 

The hybrid and thin disk models both End p = 1.0 ± 0.3 for the slope of the dependence 
of the disk size on wavelength, while the Gaussian model favors a steeper (in temperature) 
slope of p = 0.7 ± 0.3. The three estimates are mutually compatible and smaller than the 
prediction of the standard thin disk model {p = 4/3). The differences for the (presnmably) 
more physical hybrid or thin disk models are small enough to represent only a modest 
inconsistency. Experiments with broader wavelength ranges are needed to provide better 
estimates of the temperatnres exponent. 
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Fig. 1.— Magnitude differences as a function of the inverse of the observed wavelength for 
Q 2237+0305 during the six epochs of observation. Image A showed chromaticity on epoch 
HJD 2452879 (a), as previously analyzed by Mosquera et ah (2009). Here we hnd that image 
A also shows evidence of chromaticity on HJD 2454056 (d) and HJD 2454404 (f), as does 
image B on HJD 2454056 (d). The shaded regions correspond to the wavelength location 
and full width of the most prominent quasar broad emission lines. The horizontal error bars 
indicate the FWHM of the hlters. 
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P 


Fig. 2.— Joint probability density function V{rs,p) in terms of the half-light radius R 1 / 2 , 
at rest Ai = 1736 A, and the logarithmic slope p (r^ oc A^). The contours for each model 
correspond to the l-cr and 2 -ct levels for one parameter. The Gaussian, hybrid and thin disk 
models are shown by the green (dashed), blue (solid) and dotted (red) colors (lines). The 
hlled squares are the Bayesian estimates for the expected values of R 1/2 and p. 
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Fig. 3.— Joint probability density function V{rs,p) in terms of the half-light radius i?i/ 2 , 
at rest Ai = 1736 A, and the logarithmic slope p (r^ oc A^) for the hybrid model. The 
separation between consecutive contours corresponds to 0.25 a, the 1-cr and 2-cr for one 
parameter contours are heavier. The blue hlled square is our Bayesian estimate for the 
expected value of R 1/2 and p for the hybrid model (the values for the Gaussian and thin 
disk models are also plotted as green and red hlled squares respectively for comparison). 
The open (hlled) circle corresponds to the measurement by Eigenbrod et al (2008) with 
(without) a velocity prior. Straight lines correspond to the measurements by Poindexter & 
Kochanek 2010b (dashed line), Sluse et al. 2011 (dotted-dashed line) and Mosquera et al. 
2013 (continuous line), that have no estimate on p; the associated error bars correspond to 
their ±la uncertainties. The size comparisons have been made setting the mean mass of the 
stars to (M) = 0.3 Mq. 
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Table 1. Q 2237+0305 PHOTOMETRY 


Filter niB — niA rnc — tha mriD — ruA Observation date 


Str-u 

(A 

= 

3510 

A) 

1.18+0.07 

Str-v 

(A 

= 

4110 

A) 

1.15+0.06 

Str-b 

(A 

= 

4670 

A) 

1.19+0.06 

V-band 

(A 

= 

5300 

A) 

1.19+0.04 

Str-y 

(A 

= 

5470 

A) 

1.19+0.03 

Iac#28 

(A 

= 

6062 

A) 

1.14+0.10 

Ha 

(A 

= 

6567 

A) 

1.15+0.06 

Iac#29 

(A 

= 

7015 

A) 

1.16+0.06 

I-band 

(A 

= 

8130 

A) 

1.16+0.08 

Str-u 

(A 

= 

3510 

A) 

1.13+0.08 

Str-v 

(A 

= 

4110 

A) 

1.13+0.10 

Str-y 

(A 

= 

5470 

A) 

1.10+0.10 

Ha 

(A 

= 

6567 

A) 

1.14+0.05 

I-band 

(A 

= 

8130 

A) 

1.15+0.06 

Str-u 

(A 

= 

3510 

A) 

1.22+0.09 

Str-v 

(A 

= 

4110 

A) 

1.14+0.02 

Str-b 

(A 

= 

4670 

A) 

1.22+0.05 

Str-y 

(A 

= 

5470 

A) 

1.22+0.08 

Iac#28 

(A 

= 

6062 

A) 

1.30+0.03 

Ha 

(A 

= 

6567 

A) 

1.02+0.08 

Iac#29 

(A 

= 

7015 

A) 

1.17+0.07 

I-band 

(A 

= 

8130 

A) 

1.05+0.08 

Str-u 

(A 

= 

3510 

A) 

1.33+0.10 

Str-v 

(A 

= 

4110 

A) 

1.28+0.06 

Str-b 

(A 

= 

4670 

A) 

1.37+0.05 

Str-y 

(A 

= 

5470 

A) 

1.31+0.04 

Iac#28 

(A 

= 

6062 

A) 

1.28+0.10 

Ha 

(A 

= 

6567 

A) 

1.24+0.06 

Iac#29 

(A 

= 

7015 

A) 

1.25+0.08 

I-band 

(A 

= 

8130 

A) 

1.25+0.06 

Str-u 

(A 

= 

3510 

A) 

1.43+0.03 

Str-v 

(A 

= 

4110 

A) 

1.32+0.08 

Str-b 

(A 

= 

4670 

A) 

1.39+0.09 


1.56+0.11 

1.47+0.06 

Aug 

21 

2006 

1.58+0.08 

1.46+0.06 

Aug 

21 

2006 

1.59+0.07 

1.61+0.06 

Aug 

21 

2006 

1.64+0.06 

1.56+0.06 

Aug 

21 

2006 

1.59+0.05 

1.59+0.04 

Aug 

21 

2006 

1.50+0.10 

1.51+0.10 

Aug 

21 

2006 

1.51+0.08 

1.47+0.03 

Aug 

21 

2006 

1.54+0.06 

1.48+0.06 

Aug 

21 

2006 

1.59+0.15 

1.48+0.06 

Aug 

21 

2006 

1.49+0.05 

1.44+0.07 

Sep 

23 

2006 

1.55+0.10 

1.41+0.10 

Sep 

23 

2006 

1.39+0.10 

1.34+0.10 

Sep 

23 

2006 

1.47+0.08 

1.38+0.06 

Sep 

23 

2006 

1.53+0.09 

1.41+0.06 

Sep 

23 

2006 

1.74+0.08 

1.74+0.10 

Nov 

17 

2006 

1.57+0.03 

1.41+0.03 

Nov 

17 

2006 

1.74+0.09 

1.63+0.05 

Nov 

17 

2006 

1.73+0.11 

1.60+0.06 

Nov 

17 

2006 

1.59+0.07 

1.62+0.04 

Nov 

17 

2006 

1.38+0.07 

1.22+0.04 

Nov 

17 

2006 

1.64+0.12 

1.43+0.05 

Nov 

17 

2006 

1.31+0.13 

1.19+0.08 

Nov 

17 

2006 

1.84+0.09 

1.44+0.09 

Oct 

15 

2007 

1.74+0.10 

1.29+0.06 

Oct 

15 

2007 

1.83+0.05 

1.37+0.04 

Oct 

15 

2007 

1.63+0.10 

1.18+0.06 

Oct 

15 

2007 

1.60+0.09 

1.29+0.10 

Oct 

15 

2007 

1.63+0.10 

1.25+0.06 

Oct 

15 

2007 

1.63+0.08 

1.26+0.06 

Oct 

15 

2007 

1.58+0.06 

1.21+0.06 

Oct 

15 

2007 

2.06+0.09 

1.48+0.04 

Oct 

31 

2007 

1.89+0.10 

1.37+0.09 

Oct 

31 

2007 

2.02+0.10 

1.46+0.07 

Oct 

31 

2007 
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Table 1—Continued 


Filter 

ruB - ruA 

me - mA 

m-D - ruA 

Observation date 

Str-y (A = 

5470 A) 

1.37±0.08 

1.98±0.10 

1.43±0.06 

Oct 31 2007 

Iac^28 (A = 

6062 A) 

1.33±0.09 

1.97±0.11 

1.45±0.08 

Oct 31 2007 

Ha (A = 

6567 A) 

1.34±0.06 

1.82±0.10 

1.40±0.08 

Oct 31 2007 

Iac^29 (A = 

7015 A) 

1.24±0.06 

1.67±0.10 

1.21±0.02 

Oct 31 2007 

I-band (A = 

8130 A) 

1.22±0.08 

1.69±0.12 

1.22±0.08 

Oct 31 2007 






